#### Strata Table #### 
strata <- lapop12 %>% 
  group_by(ur, estratopri, estratosec) %>% 
  summarise(n=n(), 
            trial = sum(trial))
strata$prop_treat <- strata$trial/strata$n
xtable(strata)

